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ABSTRACT 

Aims. We present the results of a spectroscopic analysis on three young embedded sources (HH26 IRS, HH34 IRS and HH46 IRS) 
belonging to different star-forming regions and displaying well developed jet structures. The aim is to investigate the source accretion 
and ejection properties and their connection. 

Methods. We used VLT-ISAAC near-IR medium resolution (R ~ 9000) spectra (H and K bands) to derive, in a self-consistent way, 
parameters like the star luminosity, the accretion luminosity and the mass accretion rate. Mass ejection rates have also been estimated 
from the analysis of different emission features. 

Results. The spectra present several emission lines but no photospheric features in absorption, indicating a large veiling in both H 
and K bands. In addition to features commonly observed in jet driving sources ([Feii], Hi, Hi, CO), we detect a number of emission 
lines due to permitted atomic transitions, such as Nai and Tii that are only 2-5 times weaker than the Bry line. Some of these features 
remain unidentified. Emission from Nai 2.2yum doublet is observed along with CO(2-0) band-head emission, indicating a common 
origin in an inner gaseous disc heated by accretion. We find that accretion provides about 50% and 80% of the bolometric luminosity 
in HH26 IRS and HH34 IRS, as expected for accreting young objects. Mass accretion and loss rates spanning lO^^-IO"** Moyr"' 
have been measured. The derived MiosJM.^cc is ~0.0I for HH26 IRS and HH34 IRS, and >0.1 for HH46 IRS. These numbers are in 
the range of values predicted by MHD jet launching models and found in the most active classical T Tauri stars. 
Conclusions. Comparison with other spectroscopic studies performed on Class Is seems to indicate that Class Is actually having 
accretion-dominated luminosities are a limited number. Although the analysed sample is small, we can tentatively define some criteria 
to characterise such sources: they have A^-band veiling larger than 2 and in the majority of the cases present IR features of CO and 
Nai in emission, although these do not directly correlate with the accretion luminosity. Class Is with massive jets have high L^cc/Lboi 
ratios but not all the identified accretion-dominated objects present a jet. As suggested by the SEDs of our three objects, the accretion- 
dominated objects could be in an evolutionary transition phase between Class and I. Studies of the kind presented here but on larger 
samples of possible candidates should be performed in order to test and refine these criteria. 
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K*" 1 . Introduction 12006 *). It is nevertheless important to test accretion/ejection mod- 
els in young sources at earlier stages of evolution, when accre- 
The process of mass accretion accompanymg the formation of ^^^^ dominates the energetics of the system and thus the mech- 
H solar type stars is always associated with mass ejection in form ^^^^^ ^^^^^^^ ^„g^j^j. momentum is expected to be more ef- 
of coUimated jets, that extend from few AU up to parsecs dis- ^^^^^^ ^^^^ ^^^^ ^„ interesting sample of objects are Class 
tance from the exciting source. According to an established class j ^^^^^^^^ j ^ ^1^^^ embedded stars characterised by a 
of models (Konigl & Pudritz 2000; Casse & Ferreira 2000), ac- ^jg^pjy ^^^^^^ ^p^^j^^j ^^^^gy distribution (SED) between 2 
cretion and ejection are indeed intimately related through the ^^^^^jjy considered younger than visible T Tauri 
presence of a magnetised accretion disc: the jets carry away ^^^.^ ^^^^ ^lass II sources). However, the high extinction pertain- 
the excess angular- momentum, so that part of the disc mate- ^^^^^ ^^jects strongly limits the measurement of their stel- 
rial can move toward the star The efficiency of this process is accretion properties, needed to prove that they are indeed 
measured by the ratio between the mass ejection and mass ac- ^ ^^^^^ ^^^^^^ accretion with respect to Class II sources, 
cretion rates, and depends on the jet acceleration mechanism at g^^g^^j assumption so far applied has been that most of the 
work. Measurements of such an efficiency have been so far ob- bolometric luminosity of Class I objects is due to accretion, 
tained only for classical T Tauri stars, whose accretion proper- 
ties are rather well studied through the observation of the excess Recently, however, thanks to the use of high dispersion sen- 
continuum emission at optical and UV wavelengths Values of sitive instrumentation, it has become possible to define the stel- 
Mioss /M ace in the range 1 - 1 0% have been found by diff^erent stud- lai" properties of small samples of Class I stars through their 
ies (e.g. lHartiganet ai]|l995l: IWoitas et all 120051: iFerreira et all weak photospheric lines detected in the optical scattered fight 

(IWfi ite & Hillenbrand 2004 ) and in the near-IR direct emis- 

Send offprint requests to: Simone Antoniucci, email address: sion (Greene & Lada|2002|;|Nisini et al.|2005a|;LDoppmann et aQ 

antoniucci@oa-roma.inaf.it 2005). Such studies have shown that the characteristics of Class 

* Based on observations collected at the European Southern I objects vary in fact significantly. In particular, the accretion lu- 

Observatory, Chile (ESO Programmes 072.C-0176(A), 074.C-0235(A)) minosity may span from few percent up to 80% of the bolomet- 
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ric luminosity; these findings show that not all sources defined 
as Class I are indeed actively accreting objects and suggest that 
a c lassification based on d i ff'eren t criteria is indeed required. 

IWhite & Hillenbrandl ( |2004 . notably, derived from their 
analysis of the scattered light spectra of Taurus-Auriga sources, 
that the average M\oss/Mncc for Class I of their sample is larger 
than for Class II and close to unity. They interpret this result as 
due to an observational bias induced by the effect of disc orien- 
tation in their optical spectra; if the Class I are seen prevalently 
edge-on, then the extended region emitting the forbidden lines 
from which they estimate Mioss is seen more directly than the 
obscured stellar photosphere. It is clear that such kind of biases 
can be minimised by performing spectroscopic observations di- 
rectly in the IR, where features originating in the photosphere, in 
the accretion region and in the jet can be simultaneously detected 
by instrumentation which is sensitive enough. 

In this framework, we report here the results of near-IR spec- 
troscopic observations at medium resolution of three embedded 
sources (HH34 IRS, HH26 IRS and HH46 IRS) displaying well 
developed jet structures and having a spectral index between 2 
and 10//m, typical of Class I objects. We have derived accretion 
and ejection parameters of these sources through the analysis of 
the different features detected on the spectra. The goal is to study 
how much of their energy is due to accretion and to investigate 
the efficiency of the ejection mechanism. 

We describe the sample and the observations in Sec. 2 and re- 
port the results in Sec. 3; in Sec. 4 we present the procedures ap- 
plied for the analysis of the data and infer the physical properties 
of the objects and their jets. These results are then discussed in 
Sec. 5, where a comparison with similar sources analysed in pre- 
vious studies will be also made. Main conclusions of our work 
are summarised in Sec. 6. 



2. Description of the sample and observations 

2.1. The sample 

Our sample consists of three sources, namely HH26 IRS, HH34 
IRS and HH46 IRS, selected among sources having a positive 
spectral index between 2 and 10 fim and showing well defined 
jet structures observed in the near-IR very close to the source 
itself. The basic observational properties of the targets are listed 
in TabH 

HH34 IRS (located in the L1641 cloud, at a distance of 460 
pc) is the exciting source of the spectacular parsec-scale HH34 
flow dBallv & Devind 19941). I t was initially d etected as a fain t 
optical star bv lReipurth et al.l (Il986h . but Reip urth et all ( l2000l) 
showed that the optical emission is coming from a very compact 
reflection nebulosity, while the source itself is an embedded IR 
object displaced at about 0.1 arcsec from this nebulosity. HH34 
IR S has been dete cted by IR AS and at sub-mm /mm wavelengths 
bv lReipurth etal J (1993) and I Johnstone et all (|^006), who show 
that a dense cloud core surrounds the object. The HH34 jet 
has been studied in detail both in the optical dReipurth et al 



20021; iBacciotti & Eisloffel iT999h and in die IR jPodio eTal 



20061; iTakami et al.ii2006i; iGarcia Lopez et alj|2007l) . where the 



H2 molecular counterpart of the atomic jet has been also de- 
tected. Jet proper motion analysis has shown that the jet is ori- 
ented with a n angle ra nging fro m 20 to 30° to the p lane of the 



splane ot 
thlll992h . 



sky (Eisloffel & Mundt i 1992; He aflicote & Reipurth 

HH26 IRS (located in the LI 630 molecular cloud at a dis- 
tance of 450 pc) is an embedded object, recognised to be the 
exciting source of th e HH26A/C / D cha i n and of the associated 
molecular outflow bv lDavis etaTI (ll997l) . lDavis et al.l(l20b2i) dis- 



covered that this source is also driving a small-scale H2 jet with 
a length of a few arcsec. An IR spectrum of HH26 IR from 2 to 
4 fim has been obtained by Simon et al. (2004), showing emis- 
sion line features, from H2, HI and CO, typical of active young 
sources. Observations at s ub-mm/rnm wav elen gths of this ob- 
j ect h ave been obtained bv lLis et akl d 19991) and lJohnstone et all 

daooiD. 

HH46 IRS is located in an isolated Bok globule close to the 
Gum nebula (at a distance of 450 pc). The object actually con- 
sists of an embedded young binary system of 0.26"separation 
dReipurth et al.ll2000l) that illuminates a reflection nebula from 
which the HH46/47 system of Herbig Haro objects emerges. 
At infrared wavelengths, the redshifted counter-jet has been de- 
tected penetratin g through the refle ction nebula down to the 
exciting source dEisloffe l et al.'4994'). HH46-IRS is associated 
with the IRAS object IR AS08242-5050 and has been detected 
at millimetre by Reipu rth et al.l dl993 ) . A spectrum at mid-IR 
wavelength of HH46-IRS ha s been recently obtained by Spitzer 
dNoriega-Crespo et al.l l2004l) . Proper motion measurements of 
the HH46/47 objects indicate the HH46 IRS je t is in clined by 
-34° to the plane of the sky. (.Eisloff^el & Mundllll994 . 

2.2. Observations 

The observations were performed using the VLT-UTl ISAAC 
spectrograph in SW mode on 28, 29 and 30 December 2004. 
Three spectral regions have been investigated: one in the H band 
(1.57-1.65 jum for HH34 IRS and HH46 IRS and 1.61-1.69 /urn 
for HH26 IRS), and two contiguous segments in the K band 
(2.08-2.30 jum). These spectral regions have been chosen be- 
cause here he important emission lines tracing accretion and 
ejection (e.g. Hi and [Feii] ), and also because of the presence 
of several diagnostic photospheric absorption features which, if 
detected, are suitable for spectral classification of the forming 
star 

All spectra were obtained using a 0.3"slit, providing a final 
nominal spectral resolution of about 10000 in the H band and 
8900 in the K band. The slit has been aligned parallel to the jet. 
All the work on the raw images has been performed using the 
IRAF software package, following the standard procedures for 
bad pixel removal, flat-fielding and sky subtraction. Spectra of 
standard stars have been acquired at airmasses similar to those 
of the scientific spectra and used, after removal of any intrinsic 
line, to correct for the telluric absorptions and to calibrate the 
flux scale. 

Observations of science targets and standards were per- 
formed with similar seeing conditions (~0.6"for HH34 and 
HH46 IRS and ~l"for HH26 IRS), so that flux losses are as- 
sumed to be comparable in both cases for point-like objects and 
no correction has been applied to the spectra of the sources (for 
jets see Sec. 14.21 ). 

From the calibrated spectra we derived an estimate of the 
magnitudes in the H and K bands, which we list in Tab. [T] to- 
gether with the 2MASS magnitudes. The agreement is very good 
for HH34 IRS while for the other sources the differences are 
<0.6 mag except for the H magnitude of HH46 IRS; in this case, 
however, the 2MASS H magnitude appears to be strangely high, 
being almost equal to the K magnitude, providing aJ-H colour 
unusual for this kind of sources. 

Wavelength calibration has been performed using the atmo- 
spheric OH emission lines whenever it was possible, or, alter- 
nately, using Xenon/ Argon lamp lines, always refining the align- 
ment against known atmospheric absorption features. The de- 
scribed procedure leads to a calibration error that we estimate of 
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Table 1. The observed targets and their main observational properties. 



Source 


R. A. (2000) 


DEC.(2000) 


OTj 


niH 




D 












(mag) 


(mag) 


(mag) 


(pc) 




(Lo) 








2MASS<'> 


2MASS*'> This work<'> 


2MASS'*' This work*** 








HH 26 IRS 


05 46 03.9 


-00 14 52 


16.77 


14.07 14.6 


11.88 12.3 


450 


2.01 


4.6-9.2 


HH 34 IRS 


05 35 29.9 


-06 26 58 


15.06 


13.60 13.5 


12.38 12.4 


460 


1.14 


12.4-19.9 


HH 46 IRS 


08 25 43.9 


-51 00 36 


14.20 


12.88 14.8 


12.72 13.4 


450 


1.96 


<15.0<"' 



Notes, (a) The spectral index a = (iLog(/lF,i)/'^Log(/l) is calculated between 2 and 10 yUm. (b) The magnitude values in the H and 7 bands have 
been esti mated from the calibra ted spectra, (c) Total luminosity: the source is a binary (see text for details). References. (1) from the 2MASS 
catalogue dSkrutskie et"ai]|2006h . 



the order of 0. 1 A for the H band and the first K band segments, 
and 0.2 A for the other K band spectral region investigated. 

Velocity information on the observed lines (Vlsr and AV) 
has been obtained through a Gaussian fit. The uncertainty on the 
line widths is of about 10 km s"', as we estimate considering the 
typical difference between the fitted FWHM and the one directly 
measured on the spectra. Concerning radial velocities, the con- 
tribution Vcioud of the star-form ing clouds hos ting the sources, 
i.e. 11 km s " for HH26 IR (Gi bb et al.llT995h . 8.5 km s ' for 
HH34 I RS (jA ngladaet al. 1995) and 8 km s"' for HH46 IRS 
in Vela dChernin & Massonlll99lh . has been subtracted from the 
measured Vlsr- 



3. Results 

3.1. Observed spectral features 

The continuum-normalised spectra of the three sources are 
shown in Fig. [HIS] and [3] The spectra are characterised by the 
presence of several emission lines, while there seems to be no 
evidence of absorption features from the central object photo- 
sphere. Depending on the rms noise measured on the spectra, 
we estimate upper limits on the line/continuum ratio of the un- 
detected photospheric lines between ~ 5-30 %. These ratios give 
an indication of the minimum absorption line equivalent width 
measurable for a given spectrum. Hence, a lower limit to the 
veiling can be inferred using the relationship r = EW/EW - 1, 
connecting the veiling r to the intrinsic (EW) and measured 
(EW) equivalent width of a line. Assuming a late K spectral type 
for the sources, being typical intrinsic equivalent widths around 
1 .5A for the stro ngest abso rption features in the K band in this 
type of stars (e.g. [Wallace & Hinkle 1997), one derives veilings 
rK > 1,2,5 for HH46 IRS, HH26 IRS and HH34 IRS, respec- 
tively. 

Among the emission lines, there are several H2 and [Feii] 
transitions coming from the jet, as testified by the fact that these 
lines are not confined to the source but appear extended in the 
jet direction in the spectral images. [Feii] lines trace regions 
with moderate ionisation, temperatures between 8000 and 15000 
K and electron densities up to 10^ cm"^, which are commonly 
found in the initial part of a jet, close to the source. Rovibrational 
H2, on the other hand, is excited in a molecular gas at tempera- 
tures of the order of 2000-3000 K and total densities of -10"*- 
10^ cm-l While HH34 IRS and HH46 IRS show both [Fen] 
and H2 lines, in HH26 IRS only these latter are detected. This 
could be an indication that the jet has low excitation conditions 
or it is composed mostly of molecular material. Similar molec- 
ular jets have bee n found in a small number of Class I sources 
dPavis et alJl200lh . 



A full account of the jet physical properties derivable 
from the detected emission lines will be g iven elsewhere 
jGarcia Lopez et al]|2007l : ICaIzoletti et al.l2 007, in preparation). 

Different permitted emission lines are also detected, the most 
prominent being the H i recombination lines of the Brackett se- 
ries (see Fig. |5]l and the Na i doublet lines around 2.20 fj.m (Fig. 
ISll, which are observed in HH34 IRS and HH26 IRS. Such lines, 
whose emission is confined toward the source, trace high density 
gas excited either in the accretion region or at the base of stel- 
lar/disc winds. The H i lines have large widths as is usual in these 
kind of sources (e.g. Nisini et"ani2005 a): the Bry line in HH34 
IRS has a width of order 240 km s~' while in the other sources it 
is of order 150 kms"'. The peak of the emission is just slightly 
redshifted, but the profile is not always symmetric, as it extends 
more in the blue in HH26 IRS. In the case of HH46 IRS the 
line profile is difficult to interpret since the line is merged with 



Fig. 1. The H and K band normalised spectra of HH26 IRS. 
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Fig. 2. The H and K band normalised spectra of HH34 IRS. 



Fig. 3. The H and K band normahsed spectra of HH46 IRS. 
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another prominent emission feature whose identification is not 
certain, although the wavelength (2.165 fim) is compatible with 
He I. The velocity shift of about -200 kms ' (similar to the one 
measured on Fe ii) could indicate a Bry high velocity component 
associated with the Fe ii jet but confined on the source and not 
extended along the jet axis. Such a high velocity component is 
however not detected in the other Brackett series lines. 

Bry equivalent widths measured on the sources span from 
-4.8 to -8.9 A(seeTab.[8ll. 

Band-head emission from CO is also detected in the sources, 
even though the emission is much weaker in HH46 IRS; this is a 
tracer for dense molecular gas at T ~2000-3000 K, usually inter- 
preted as a sign of the presence of either an active accretion disc 
or of a neutral wind (e.g. Carr 1989; Carr et al. 1993; Naiita et al. 
119961) . 

Finally, we detect other lines for which we provide tentative 
identifications, namely Tii lines at 2.179 and 2.190 //m (visible 
in HH26 IRS and also in HH34 IRS) and other features that re- 
main unidentified (e.g. at 2.239 //m in HH26 IRS and at 2.240, 
2.281 ^lm in HH34 IRS). The Tii lines (that are well detected 
especially HH26 IRS) are broader than H2 lines, while they ba- 
sically have the same width as the Na i doublet lines (around 80 
km s"'), thus suggesting that they probably originate in the same 
region around the central source. 

A list of all the observed emission features having S/N 
greater than 3 and line-widths at least equal to the instrumen- 
tal resolution is reported in Tabs [3] |4] and |5] together with 



their fluxes and tentative identifications based on the and 
Atomic Line LisQ 

Radial velocity information for the brightest lines of the de- 
tected series (H i Bry, H2 at 2. 12 yum, Fe 11 at 1 .64 jum and Na i at 
2.2 ^m) is given in Tab.|6] 

3.2. Bolometric luminosity 

The total luminosity Lboi of the sources has been derived by in- 
tegrating the observed spectral energy distributions that include 
measurements spanning from the near IR to the mm region (see 
table |2]i. The calculation is performed starting from the J band 
value and considering straight lines (in the Log(/l)-Log(/lF^) 
plan) between available SED points; a final correction at the 
longest wavelengths is applied assuming that the emission 
is decreasing as l/A^ after the last available observation at 1.3 
mm. The derived value might be subject to a bias because of 
the possible presence of unseen companions contributing to the 
observed flux or contamination from other stars in the field, es- 
pecially for the IRAS measurements at 60 and 100 microns. For 
this reason, we considered recent Spitzer-MIPS images to ver- 
ify the crowding of the fields at mid IR wavelengths and derive 
more accurat e flux estimates. Results for HH46 IRS at 24 jim 
are given by dNoriega-Crespo et al. 2004. while for HH26 IRS 
and HH34 IRS we have used available archive images at 24 and 
70 jum. Spitzer photometry at 24 and 70 microns was considered 
instead of the respective IRAS fluxes at 25 and 60 microns. 



http://physics.nist.gov/PhysRefData/ASDl 
httpV/www.pa.ulcy.edu/ peter/atomic databases 
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Fig. 4. Spectral Energy Distributions of the three sources, derived from 
the values reported in Tab. [2] Points from different instruments are dis- 
played with different symbols (2MASS: crosses, MSX: triangles, IRAS: 
squares, SPITZER: diamonds, far IR/mm measurements: stars. 
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For HH26 IRS, Spitzer photometry shows that contamina- 
tion from nearby sources is indeed very likely in the IRAS mea- 
surements. In fact, the calculated 70 ;um flux is much lower than 
the IRAS 60 jim point. A luminosity of 9.2 Lq was derived con- 
sidering the IRAS 100 jim point that however is likely to be over- 
estimated. Neglecting the 100 point we obtains 4.6 Lq. These 
results are to be regarded as an upper and lower limit to the total 
luminosity of the source. 

HH34 IRS does not appear in the IRAS point source cat- 
alogue, probably owing to source confusion problems ; how - 
ever, IRAS fluxes were extracted by Cohen & Schwartz, (|1987|) . 
Actually, the Spitzer images show that important contamination 
from other sources is not likely, as the source appears to be quite 
isolated in the field. On the other hand, the SED displays a very 
rapid increase at lOOyum with respect to the 60/im point. The 
derived Spitzer fluxes are lower than the IRAS fluxes, although 
the discrepancy is not as big as in the case of HH26 IRS. The 
derived total luminosity is about 19.9 Lq considering the IRAS 
points and 12. 4 Lq neglecting the 100 fim measurement. 

[Reipurth et^lj(|2000) found HH46 IRS to be a binary using 
HST/NICMOS observations. The binary flux ratio in the band of 



Table 2. Photometric measurements available in the literature. 



Measurements 


rlrlzu IKo 


UTJ1A TDQ 


III 1 ,1 /: TD C 

rlrl4D IKo 


zMAho J (mag) 


16.77 


15.06 


14.20 


2MASS H (mag) 


14.07 


13.60 


12.88 


2MAoo K (mag) 


1 1 CO 

1 1.00 


IZ.JO 


1 O TO 

IZ. /Z 


iKAb Iz ^m (Jy) 


1.93 


0.65 


0.81 


iIv/\o ZJ flLll yjy ) 


J. I J 




U.J 1 


Du /Jill \ jy ) 


ZU.O / 




ZD. J 1 


TT? A *v 1 no (#m 
iK/\o lUU fim yjy) 


D / .yj 




JO.Z / 


iVloA o.Zo flLli {jy) 








iVloA i'4.0j flUl yJy) 








iVlozV Z. 1 . JU ^111 yjy ) 




J. JH- 




SPITZER 24 yum (Jy) 


l.S***' 


6.7*«> 


3.9*" 


SPITZER 70 yum (Jy) 


7.9(8) 


23.3<**' 




350 yum (Jy) 


6.30^^' 






800 fim (Jy) 






l'.12<*' 


850 fim (Jy) 


l'.19<2) 


3'.4(" 




1 . 1 mm ( Jy) 






O.34W 


1.3 mm (Jy) 


O.32P' 


0,42(6) 


0.1 8<*' 


3.6 cm (mjy) 




0.16« 





Re ferences. ( 1) Johnstone & Bally (2006); (2) Johnstone et al. (2001^ 
(3)lLis et al.llfr999): (4) Correia et al, (1998); (5).Rodriguez & Reipurtl] 
{WM) : (6) lReipurthetai] ( IT993h ; (7) lNoriega-Crespo et al.l ( 120041) : (8) 
Spitzer archive. 



observation is about 0.7, but it is impossible to determine which 
of the two components is responsi ble for the jet. Also in this 
case, the Spitzer flux at 24 microns ( iNoriega-Crespo et alj|2004l) 
is lower than the one measured by IRAS. The inferred total lu- 
minosity of the binary is 15.0 L0, which of course is an upper 
limit to the bolometric luminosity of the HH46 exciting source. 



Fig. 5. H I line profiles for (from left to rigth) HH26 IRS, HH34 
IRS and HH46 IRS. 



HH26 IR HH34 IRS HH46 IRS 




Vls» - V.„„, (km s-') V,s, - V„„„, (km s"') V,s, - V„„„, (km s"') 
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Fig. 6. Nai line profiles for (from left to rigth) HH26 IRS and 
HH34 IRS. 

HH26 IR HH34 IRS 



Na I _ X ,.2- Na 




Vlsr - V„„„ (km s-') V,s„ - V,„„, (km s"') 



Table 3. Emission features in HH26 IRS. 



A F±AF 
(pm) (lO^^ergs"' cm"2) 


FWHM 

(kms-') 


ID 


1.6116 1.2 ±0.7 


140 


Br 13 


1.6413 3.6 ±1.6 


235 


Br 12 


1.6530 1.0 ±0.3 


100 


? 


1.6756 0.7 ± 0.2 


45 


7 


1.6812 2.2 ±0.5 


100 


Br 11 


1.6877 1.9 ±0.3 


60 


H2 1-0 S(9) 


2.0874 2.4 ± 0.5 


70 


7 


2.1216 69.0 ±0.4 


50 


H2 1-0 S(l) 


2.1541 6.9 ±1.2 


210 


H2 2-1 S(2) 


2.1663 26.0 ±0.8 


140 


Bry 


2.1791 4.6 ±0.5 


80 


Til 


2.1906 3.4 ±0.6 


80 


Til 


2.1905 3.1 ±0.5 


70 


7 


2.2012 2.3 ± 0.4 


60 


H2 3-2 S(3) 


2.2065 4.5 ± 0.6 


80 


Nai 


2.2092 7.8 ± 0.6 


95 


Nai 


2.2231 15.9 ±0.4 


60 


H2 1-0 S(0) 


2.2391 2.3 ± 0.5 


65 


7 


2.2475 6.6 ± 0.4 


75 


H2 2-1 S(l) 


2.2940 200.2 ± 1.5 




CO (2-0)<a) 



Notes. Fluxes are not corrected for extinction, (a) measured up to 2.305 
fim. 



3.3. Source classification 

The near-IR photometry and the availability of Spitzer data al- 
low us to refine the classification of the three sources based on 
their SEDs. The derived ^^-[24] colours are 10.3 mag for HH26 
IRS, 12.5 mag for HH34 IRS and 13.3 mag for HH46 IRS. These 
values indicate that the three sou rces can be indeed c onsidered 
as Class I or younger, following iRebuU et all (l2007l in press) 
who classify as Class I sources those having a K-Q colour larger 
than 8.3. The near-IR J - H and H - K colours of HH34 IRS 
and HH26 IRS are instead typical of embedded T Tauri stars 
(|]VIeyeretal.lll997h . The HH46 IRS colours are rather peculiar, 
since the source falls on the left of the reddened main-sequence 
star locus. As we mention in Sect. 14.1! however, the displace- 
ment of the sources in this colour-colour diagram is not indica- 
tive of the true amount of extinction, since the presence of a large 
quantity of scattered light has the efi'ect of making the star appear 
bluer and less reddened (see e.g iMassi et al.l [19991 IStark et aTl 
|2006|) . Contamination due to light scattered by the circumstellar 
disc is indeed very likely in all our sources, that have estimated 
disc inclination angles, with respect to the line of sight, of 20- 
30°. 

Indication that the sources are in an early stage of evolution 
is al so given by their relatively massive dust envelopes (0.2-1 
M0 iReipurth et al.lll993l: iLis et al.lll99"9l) . Indeed, the ratio be- 
tween Lboi and the sub-mm luminosity measured longward of 



Table 4. Emission features in HH34 IRS. 



A 


F±AF 


FWHM 


ID 


(jum) 


(10""^ ergs"' cm-^) 


(kms-') 




1.5773 


1.5 ±0.3 


85 


7 


1.5888 


4.2 ± 0.7 


265 


Brl4 


1.5895 


2.3 ± 0.3 


90 


7 


1.5996 


37.6 ± 0.4 


100 


[Fell] 


1.6115 


4.3 ± 0.6 


180 


Brl3 


1.6413 


7.9 ±0.1 


290 


Brl2 


1.6437 


156.0 ± 1.0 


100 


[Fell] 


2.1220 


22.0 ±0.1 


15 


H2 1-0 S(l) 


2.1329 


8.7 ± 0.5 


80 


[Fell] 


2.1545 


1.1 ±0.3 


60 


H2 2-1 S(2) 


2.1663 


39.0 ± 0.1 


240 


Bry 


2.1793 


2.5 ± 0.5 


55 


Til 


2.1902 


3.5 ± 1.1 


170 


Til? 


2.2016 


1.2 ±0.3 


45 


H2 3-2s(3) 


2.2065 


8.3 ± 0.7 


115 


Nai 


2.2092 


7.9 ± 0.7 


110 


Nai 


2.2234 


9.1 ±0.2 


10 


H2 1-0 S(0) 


2.2239 


34.0 ± 0.7 


95 


[Fell] 


2.2403 


2.3 ± 0.4 


55 


7 


2.2438 


6.3 ± 0.6 


85 


[Fell] 


2.2478 


2.2 ± 0.2 


10 


H2 2-1 S(l) 


2.2536 


8.8 ±0.5 


85 


[Fell] 


2.2807 


2.1 ±0.6 


75 


7 


2.2940 


184.0 ± 1.5 




CO (2-0)<a) 



Notes. Fluxes are not corrected for extinction, (a) measured up to 2.305 
fim. 

Table 5. Emission features in HH46 IRS. 



A 


F±AF 


FWHM 


ID 


(fim) 


(10""^ ergs"' cm-^) 


(kms-') 




1.5889 


5.1 ±0.8 


310 


Br 14 


1.5989 


18.3 ±0.3 


60 


[Fell] 


1.6114 


2.7 ± 0.4 


150 


Br 13 


1.6387 


12.1 ±0.5 


100 


7 


1.6411 


21.1 ±0.5 


120 


Br 12 


1.6429 


173.2 ± 0.4 


65 


[Fell] 


2.0847 


1.3 ±0.4 


80 


7 


2.1218 


14.8 ±0.3 


45 


H2 1-0 S(l) 


2.1321 


2.4 ± 0.4 


110 


[Fell] 


2.1647 


3.0 ±0.3 


80 


7 


2.1662 


7.4 ± 0.6 


175 


Bry 


2.2231 


20.3 ± 0.7 


80 


H2 1-0 S(0) 


2.2432 


5.1 ±2.0 


350 


[Fell] 


2.2477 


2.5 ± 0.5 


50 


H2 2-1 S(l) 


2.2526 


3.4 ± 0.9 


85 


[Fell] 


2.2936 


44.8 ± 0.9 




CO (2-0)<a) 



Notes. Fluxes are not corrected for extinction, (a) measured up to 2.305 
fim. 



Table 6. Radial velocities ( Vlsr- Vcioud) measured on the bright- 
est lines of the detected series. 



Source 


Hi Bry 


Fen (1.64;um) 


H2 (2.12;um) 


Nai 




(kms-') 


(kms-') 


(kms-') 


(kms-') 


HH 26 IRS 


-3 




-56 


-10 


HH 34 IRS 


-3 


-88<°> 


-11 


3 


HH 46 IRS 


-3 


-200 


-11 





Notes. We estimate an uncertainty of 5 kms-' for the values liere 
reported, (a) the line is actually double-peaked, with a blue peak 
Vlsr - Vdoud=-93.5 kms-' and a red peak Vlsr - Vdoi,d=-62.5 kms-'. 
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350 fim (Lsubmm), estimated from the photometric data, is rather 
small in all the sources ( -60 for HH26 IRS, ~ 70 for HH34 IRS 
and ~160 for HH46 IRS), which would indicate that they actu- 
all y satisfy the Class criterium Lboi / i-submm < 200 as defined 
bv lAndreetai] (119931) . Moreover, co mparing the observed S EDs 
with the grid of YSO S EP models bv I Whitney et alj (l2003h and 
lRobitailleetan (l200g). we derive that models for evolved Class 
sources provide better fits to the observational data than those 
for Class I/IIs, although some parameters, such as the depth 
of the silicate feature, largely deviate with respect to the ob- 
served values. In any case, these findings strongly suggest that 
the sources are indeed very young and are not more evolved 
(Class II) objects. 



4. Analysis and physical parameters 

4.1. Extinction, Accretion and steiiar parameters 

The fraction of the total luminosity due to accretion (L^cc) must 
be d etermined in order to get an estimate of the mass accretion 
rate. iNisini et"aLl ( 12005 a l) derived this quantity for a sample of 
young sources directly determining the spectral type and veil- 
ing from the analysis of the absorption features from the stel- 
lar photosphere, then inferring the stellar parameters and lumi- 
nosity L«, from which the accretion luminosity can be obtained 
(Lace = Lboi -L»). However, our spectra show no evidence of well 
detected (S/N > 3) absorption features from the photosphere, in- 
dicating the presence of a strong veiling, so that an alternative 
method has to be used to derive the accretion luminosity. 

We have th e refore considered the relationship found by 
iMuzerolle et al.l (11998a ') on the basis of observations of T 
Tauri stars, in which Lace is directly related to the (extinction- 
corrected) H I Bry flux: 

Log ^ =(1.26) Log ^+(4.43) (1) 
L© Lq 

Such an approach implies considering this relation still valid for 
the more embedded sources of our sample; the only observa- 
tional evidence for such an assumption has been given so far 
for thr ee sources in the R CrA star-forming region jNisini et al.l 
l2005ah . 

In order to retrieve absolute Bry luminosities from observed 
fluxes, it is necessary to correct them for the extinction and thus 
to measure the value of Ak toward the sources. 

Estimates of the extinction can be inferred adopting standard 
methods based on the analysis of the spectral features. For HH46 
IRS, considering the silicate absorpt ion feauture at 9.7 fim oh- 
served by Spit zer dBoogert et al.ll2()04) and using the relation by 
lMathisl(ll998h . we obtain an value betwe en 32 and 38 mag, 
corres ponding to an 3.5-4.2 adopting the iRieke & Lebofskvl 
(11985') extinction law. In the case of HH26 IRS, measurements 
of the ice absor ption feature a t 3 yurn (Siinon et al. 2004) and the 
relationship bv iTanaka et al.l d 19901) provide an Ay value of the 
order of 29 mag (A^ ~ 3.2). 

We note, however, that the ISM relationships between dust 
features and extinction, derived for the ISM, may not be strictly 
valid for dust in the circumstellar environment of YSOs (e.g. 
IChiar et al ] l2007l) . 

Finally, for HH34 IRS only a lower limit of Ay =7 mag can 
be tak en adopting the value recently estimated by Podio et aO 
( |2006|) from the ratio of [Fe ii] lines excited in the jet region at a 
distance less than ~l"from the source. 
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Given these problems, since a reliable determination of the 
extinction is crucial in order to derive the accretion luminos- 
ity from Bry, we have tried to use a different and independent 
method to compute "self-consistent" Ak, based on the observed 
magnitudes and on reasonable assumptions on the stellar proper- 
ties of the objects. Indeed, the intrinsic stellar luminosity is given 
by: 

LogLJLQ = -0.4(Mboi - Mb„,,o) (2) 
where 

Mboi = BC + Mk + {V - K), (3) 

The absolute magnitude Mk in Eq. |3] depends on the observed 
magnitude mg, the extinction A/j-, the veiling and the distance 
modulus DM: 

Mk = mk + 2.5Log(l +rK)-AK- DM (4) 

while the values of the bolometric correction BC and the intrin- 
sic stellar colour {V - K)^ depend on the spectral type and age 
of the objects. Thus, if we assume that our sources are late-type 
stars located on the birthline and use the parameters provided by 
a set of evolutionary models, we can obtain an estimate of Ak 
from the observed ^T-band magnitude and veiling of the objects. 

On this basis, we have adopted a procedure allowing us to 
derive the different physical quantities in a self-consistent way. 
Namely, we have varied the Ak until we have found a value able 
to reconcile the stellar luminosity given by Eq. 2 with the value 
given by Lboi - Lace, where Lace is derived from the extinction- 
corrected (by the considered Ak) Bry line emission using Eq. 
1. 

For example, we can easily derive that an extinction value 
Ay -29 mag for HH26 IRS (such as the one derived from the 
ice feature) would lead to a stellar luminosity which is smaller 
than the value given by Lboi- Lace (A L, ~ 2.1Lq), pointing to 
an inconsistency in the independently inferred A^ and L, values. 

We have carried out the described pro cedure assuming the 
birthline location of Pa lla & Stahled d 19931) and considering the 
colours and bolometric correction given by [Siess et al. (200^ 
to convert L, into Mk (Eq. 2 and 3); thus, we have found the 
extinction providing "consistent" L, values. 

We derive for HH26 IRS an Ak =4.1 mag with a spectral 
type K7 (L, = 3.7 Lq, /?,=4.1 Rq, M»=0.6 Mq) and for HH34 
IRS an Ak - 4.9 mag with a spectral type MO (L, = 2.9 Lq, 
/?»=3.9 Rq, M.=0.5 Mq). 

In the case of HH46 IRS we know that the source is a bi- 
nary, but the contribution from the companion remains unknown. 
Assuming that the bolometric luminosity is equally distributed 
between the two sources, the best results are obtained consid- 
ering a K5 star with Ak - 4.9 mag (L,=6.0 L©, -4.2 Rq, 
M» = L2 Mq). If we assume that the contribution to the bolomet- 
ric luminosity of our considered target is larger with respect to 
the companion, we find that it is not possible to obtain consis- 
tent results for L, smaller than 10 Lq; exploring the possibil- 
ity of a more massive star, we can still find consistency with an 
Ak -5.3 mag and a spectral type G9/K0 (L» = 12.5 Lq, /?*=4.3 
R0, M,=2.8 Mq). However in this case the bolometric contri- 
bution of the companion (considering a derived Lace ~ 2.4 Lq) 
would be basically zero, which appears unlikely. 

In Tab. [7] we report the parameters derived for the three 
sources adopting the analysis just described. 

Providing self-consistent values of the parameters, we are 
confident that these are better determinations of the extinction 
with respect to the previous ones, also taking into account the 
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already mentioned problems that may arise using the relation- 
ships established for the ISM for dust around YSOs. 

Moreover, adopting the extinction values from spectral fea- 
tures and jet line ratios, we would infer accretion luminosities 
~1 Lq for HH26 IRS, -0.1 Lq for HH34 IRS and 0.3-0.7 Lq 
for HH46 IRS, that represent only a small fraction of the total 
luminosity of the sources (~l-20%), thus implying that most of 
the source luminosity should be contributed by the stellar pho- 
tosphere. However, this is in contrast with the early evolutionary 
state of the objects and with the large K band veiling values we 
find. 

Conversely, the accretion luminosities that we infer adopting 
our "self-consistent" extinction values are 3.2 L0 for HH26 IRS, 
13.3 Lq for HH34 IRS and 1.5 Lq for HH46 IRS. The estimated 
accretion to total luminosity ratios are therefore ~0.5 for HH26 
IRS, -0.8 for HH34 IRS, and -0.2 for HH46 IRS (see Tab|2ll, 
where we have considered a mean value of the Lboi ranges given 
in Tab.[I]for HH26 IRS and HH34 IRS. 

From these accretion luminosities we have then derived an 
estimate of th e mass accretion rate M ace, using the formula for 
disc accretion (iGullbring et al.lll998l) : 



retrieve an estimate of the total flux in the considered jet section. 
We therefore have: 



Ma, 



(5) 



where and M, are the stellar radius and mass, respectively, 
and Ri is the inner radius of the accretion disc. We consider the 
masses and radii derived from Siess models during the computa- 
tion of the extinction described abov e and assume a value of Ri 
= 5 R, for the disc inner r adius (e.g. Gullbrin get ai]|l998l) . The 
inferred mass accretion rates are listed in Tab. |7]and are of the 
order of 10"^ (HH26 IRS and HH46 IRS) and 10"^ (HH34 IRS) 
Mgyr-i. 

All previous calculations were carried out assuming the mea- 
sured lower limits on veiling. Considering higher rx values and 
reperforming the computations leads to increased Ak and Lace, 
lower L» and M,, and consequently to greater Lacc/i'boi and Mace. 
In particular, in the case of HH26 IRS, taking into account a veil- 
ing r/c = 4 we get a higher Ak ~ 4.4 with a spectral type Ml 
and Laee ~ 4.IL0; we would thus obtain Laee/iboi ~ 0.6 and 



Mace -1.6-10-" Moyr-^ 



4.2. Jet mass flux 



The jet mass flux can be derived from emission line measure- 
ments, providing that the jet is resolved and its velocity is known 
(see e.g. Hartigan et al. 1995; Nisini et al. 2005b). For our pur- 
pose of directly comparing the mass loss with the mass accretion 
rate, we need to determine the mass flux relative to jet compo- 
nent as close as possible to the central source. Jet knots far from 
the protostars may indeed be related to older episodes of mass 
ejection, which could be associated to a diff'erent mass flux. 

For HH46, we have estimated the mass flux from the lumi- 
nosity of the [Fell] 1 .64 yum line, measured on the first 3 arcsec of 
the jet length. There are no measurements of the jet diameter in 
this inner jet section, which is not optically visible: at a distance 
of -5 arcsec, a jet diameter of - 1 arcsec has been measured 
by lEisloffel & Mundtl (Il994l) . Since j et diameters increase with 
distance from the driving source (e.g. lDoupados et al.l2000l) . the 
diameter in the portion of the jet we are considering should be 
- 0.6 arcsec. Thus, in addition to the flux losses due to the see- 
ing, we have also taken into account a correction factor of 2 to 



Mioss = fimHX(nHV)xv,/l, 



and 



[Fe] 
[H] 



(6) 



(7) 



where A,, / are the radiative rate and fractional population of 
the upper level of the considered transition, [Fe]/[//] is the to- 
tal Fe abundance with respect to hydrogen (assumed 2.82 -lO"^, 
Asplund 2005) and v, and I, are the velocity and length of the 
knot, projected perpendicularly to the line of sight. We assume 
here that all Fe is ionized, a s it is expect ed in the excitation con- 
ditions of HH objects (Nisini et al.ll2002l) . The fractional popula- 
tion has been computed adopting a NLTE statistical equilibrium 
code CNisini et al.l 12002,) with an electron density «(.=3.7 10^ 
cm""^, estimated from the [Fen] 1.64/1.60 observed line ratio. 
The electron temperature has been assumed 10^ K. The tan- 
gential ve locity has been i nferr ed through the proper motion 
study of Eisloff^el & Mundtl (1 1994) . and it is equal to 170 kms 
The 1.64 line luminosity has been calculated from the ob- 
served flux, which has been dereddened assuming that the ex- 
tinction in the considered jet knot is A^ = 6.6 mag. This value 
has been estimated from the [Fen] 1.64/1.25 /vm line ratio mea- 
sured by Fernandes (2000), adopting the expected theoretical ra- 
tio of lOuinet et al.,(.1996i) . The resulting mass flux value is about 
310"^ M0 yr~'. Such a valu e is about a factor of 3 and 1 lower 
tha n the values estimated by Bacciotti & EislofFel (1999) and by 
Hartigan et al.i (1994) , respectively, through the measurement of 
the total gas density, and thus the total mass in the jet, follow- 
ing the procedure known as the BE technique (from Bacciotti & 
Eisloff'el, 1999), or adopting a shock model. It has been shown 
that the mass flux thus inferred may be overestimated due to 
the assumption that the considered jet area is filled with ga s at 
the derived density (iHartigan et al.lll994l:lNisini et al.ll2005al) . If 
the shock is caused by jet time-variability, the lack of correc- 
tion for shock compression will overestimate the time-averaged 
Mioss by a factor of about 5 dHartigan et al.lll994)) . On the other 
hand, iHartigan et al.l (1 19941) also compute a mass flux value of 
4.2- 10"^ M0yr~' from the luminosity of the [Oi] 6300A line, 
adopting a relationship similar to that given in Eq.|6] although on 
a jet region much more distant (-10-20"from the source) than 
the one we consider 

Part of the above discrepancies may arise from our assump- 
tion of solar abundance of gas-phase Fe. This assumption may 
not be valid if part of iron is still locked in grains along the 
jet. This is i ndeed what has bee n fo und in other jets (HHl 
and HH34) bv lNisini et all (l2005bl) and lPodio et al] (|2006|) . who 
showed that in the inner jet regions only a fraction between 30 
and 70% of all iron may be in gaseous form. If the same applies 
also to the HH46 jet, than our data give a mass flux up to -2-10"^ 
M0yr-i. 

The mass flux in the HH34 inner jet has been mea- 
sured by several authors. The most recent determinations were 
given by Davis et al. (2 003) , wh o used the luminosity of the 
[Fe 11] 1 .64yum line, by Podio et al.l (l2006l) . who compared diff'er- 
en t techniques employing bo th optical and IR lines, and finally 
bv lGarcia Lopez et al.l (l2007l) . who used the [Fe 11] 1 .64yum lumi- 
nosity from the same set of data we are here analysing. All these 
determinations agree on a value in the range 4-7-10"** M0yr"^ 
If we correct the determinations based on [Fe 11] line luminosity 
by the 30% of iron dust depletion found in the inner jet knot by 
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iPodio et al] ( l2006l) . we get amass flux up to -1.2 10"^ M0yr 
which is in agreement also with the determinations by Hartigan 
etal.(1994). 

The jet in HH26 is seen close to the star only in H2, indi- 
cating that the gas is mainly molecular and that the atomic com- 
ponent does not give a significant contribution to the transport 
of material outwards. We have therefore measured the mass flux 
from the luminosity of the observed 2.12//m H2 line. More ex- 
plicitly, Eq. 6 can be here written as: 



Mioss = 2// m// X iNiH2)A) x v,/l, 



(8) 



where N{H2) is the column density of molecular hydrogen and 
A is the considered emission area. Since H2 is optically thin 
and quickly thermalised, the total column density can be directly 
measured from the dereddened line luminosity: 



N{H2)A^L2.nt.n, (hvAifi)- 



(9) 



with fi given by the LTE Boltzman population at the gas ki- 
netic temperature. A temperature of ~ 2000 K can be estimated 
by the ratio of the various observed lines. We have consid- 
ered only the first knot of the HH26 microjet, which extends 
within about 2 arcsec from the central source. As in the case 
of HH46, we have applied a correction of a factor of 3 to take 
into acco unt the width of the jet a t the knot position, which 
is of ~l"( IChrvsostomou et al.l 120071) . Th e adopted line flux re- 
sults consistent with the one measured bv' Giannini et alj (|2004|) 
with a slit of 1". The tangential velocity has been estimated to 
be around 130 kms ' from the measured radial velocity of -60 
kms assuming an inclination angle of 65° with respect to the 
line of sight (Davis et al. 1997) . Here, a critical parameter is the 
reddening towards the iet. lGiannini et al.l ( l2004l) measured an Ay 
=2 on the external HH26A knot; it is however likely that the ex- 
tinction has a sharp increase close to the object, since a value of 
Ay ~38 has been estimated on-source. A method to have an esti- 
mate of Ay in the H2 emitting region is to construct a Boltzman 
diagram with the H2 lines detected on source (Table 3) and take 
the Ay value that gives the better alignment of the different tran- 
sitions on a straight line (see e.g. Nisini et al. 2002). Given that 
the considered lines do not cover a wide range of wavelengths, 
the method is not very sensitive to small Ay variations: we find 
that the best alignment is obtained with Ay values between 10 
and 20 mag. With this range of extinctions, we get H2 column 
densities ranging between 5 and 13-10"^ cm"^, corresponding to 
a mass flux ranging from ~2-10"'* to 5.3-10"'* Mq yr"'. 

The (one-sided) estimates of Mioss we have thus derived are 
summarised in Tab.[8l 



5. Discussion 

It is interesting to analyse the derived accretion and ejection 
properties and their mutual connection in the light of the results 
provided by some works that have recently investigated the char- 
acteristics of embedded protostellar sources. In particular, we 
win c ompare the results h ere presented to observations carried 
out bv^Nisini et al.l ( l2005ah (hereafter N05) on a small sample of 
Class 1 sources in the R CrA star-forming region, using the same 
instrument and observational setup. The most interesting sources 
(namely HHIOO-IR and IRS 2) displayed strong Hi emission 
lines along with well detected absorption features from which 
stellar properties and veiling could be derived. Unlike the objects 
analysed in this paper, all the CrA sources showed no signature 
of jet emission in the region closely surrounding the central star; 
only HHIOO-IR has been suggested to be the driving source of 



the Herbig Haro objects HHlOl and HH99 dHartigan & GrahamI 
1987) and of a bipolar molecular outflow (Anglada et aljl 1989b . 
which however are observed far from the central object. 

To help the comparison, we added in Tab. l7l8l also the param- 
eters of the three CrA sources found by N05 to have H i emission. 

We will also refer to two recent surveys of embedded sources 
that have investigated the nature of Class 1 sources with respect 
to the Class 11 objects, in order to verify the paradigm of the 
star formation according to which Class I are actively accretin 
objects less evolved than Class II: IWhite & Hillenbrand! (|200' 
(hereafter W04) analysed high resolution optical spectra of a 
large sample of environmentally young stars (i.e. Class I and 
Class II objects powering Herbig-Haro flows) in Taurus-Auriga. 
However, since the Class Is observed by the authors are optically 
visible, their sample is probably biassed tow ard more evolved 
objects than standard Class I sources. Doppm ann et al.l ( l2005l) 
(hereafter D05) studied high resolution near IR spectra of a large 
set of Class I and flat-spectrum sources in various nearby star- 
forming clouds and derived that Class I show in average higher 
veilings and angular velocities than what is observed in Class II 
sources, as expected for less evolved and accreting objects. 

5.1. Spectral features 

The medium resolution spectra of our sources show important 
different characteristics. In fact, both the detected emission fea- 
tures and their relative strength vary depending on the objects, 
as can be seen in Fig.[T]|2]and[3] H i, CO, and H2 emission is de- 
tected in each of three sources, which in general should be sug- 
gestive of the presence of accretion/ejection flows and of an ac- 
cretion disc. These features are usually, but not always, observed 
in embedded sources; in particular, D05 detected H i, CO and H2 
emission in 65%, 15% and 44% of the objects of their sample, 
respectively, remarking as the presence of the CO is always ac- 
companied by H I emission. This is a behaviour we find in our 
sources as well, as in those of N05. In contrast, the presence of 
permitted lines does not appear to be strictly connected with the 
simultaneous presence of jet lines from the region closely sur- 
rounding the central source. For instance, the Nai is detected 
both in HH26 IRS and HH34 IRS, but the [Fe 11] lines are visible 
only in the second source. This lack of correlation is extreme in 
the case of the CrA objects in N05, with strong permitted emis- 
sion lines and no jet signatures at all. Therefore, ejection and 
accretion signatures that are detected in the spectra do not seem 
to be strictly related. 

All these findings must denote different physical conditions 
in the close surroundings of the sources: such an indication is 
also given by the analysis of the flux ratio between Bry and Na i, 
which is not constant in the objects (see Tab. [8]). We note that 
this ratio is smaller in the jet sources with respect to the CrA 
objects; since H and Na have very different ionisation potentials 
(13.6 and 5.1 eV, respectively), a variation of the flux ratio must 
reflect a different ionisation degree in the region traced by the 
Bry and Na i. 

Na I lines are in most cases observed in Class I sources where 
CO emission is detected: 5 objects out of 8 in D05 and 2 out of 
2 in N05. On the contrary, in the majority of T Tauri stars of 
D05, both CO and Nai are observed in absorption. Noticeably, 
in our sample the CO(2-0)/Nai ratio (Tab. [8]) follows a trend 
opposite to the one of Bry/Nai ratio, i.e. it is higher in the jet 
sources. CO overtone emission is usually interpreted as coming 
from the inner gaseous disc, heated at temperatures of 2000 K 
or more by viscous accretion or magnetohydrodynamic waves 
(e.g ICarr et alJ[T993h . Nai could then be originated in an even 
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Table 7. Derived parameters for the observed targets. 



Source 


Ak 


Lace 


L, 


M. 


ST 


L-dcc I ^bol 


Mace 








(mag) 


(Lq) 


(Lq) 


(Mo) 






(10-' Mgyr-') 


(10-' Moyr-') 




HH26 IRS 


-4.1 


3.2 


3.7 


0.6 


K7 


-0.5 


8.5 


0.2-0.5 


0.02-0.06 


HH34 IRS 


-4.9 


13.3 


2.9 


0.5 


MO 


-0.8 


41.1 


0.4-1.2<''> 


0.01-0.03 


HH46 IRS 


-4.9 


1.5 


6.0 


1.2 


K5 


-0.2 


2.2 


0.3-2.0*'" 


0.14-0.90 


HHIOO-IR 


3.3+0.3 


12±2 


3.1±0.9 


0.3-0.75 


K5-M0 


0.80 


10-20 






IRS2 


2.4±0.3 


7.7+2.5 


4.3+1.5 


1.1-1.8 


K2 


0.60 


2-3 






IRS5a 


5.0±0.6 


-0.4 


1.6±0.5 


0.4-0.9 


K5-K7 


0.15 


0.2-0.3 







Notes. Columns show in ord er: K band e xtinction, accretion luminosity, stellar luminosity (assuming the star on the birthline), estimated stellar 
mass and spectral type (using fsTess et al.l (|2000) models), accretion to bolometric luminosity ratio, mass accretion rate, (one-sided) mass loss rate, 
mass loss rate to mass accretion rate ratio. The stellar luminosity and the fraction of the luminosit y due to accretion a re computed considering the 
mean value of L(,oi given in Tab.[T] The same results for three Class I sourc es of the R CrA re gion dNisini et alj2005ah are reported for comparison 
(see text for details), (a) Lower value from optical line determinations dPodio et al.l 1200 6^: upper value from [Feii] 1.64 yum line luminosity, 
assuming Fe gas abundance from lPodio et alj 1I2OO6). (b) Upper value determined assuming that only -30% of iron is in gaseous form. 

Table 8. Derived parameters for the observed targets. 



Source 


rK 


BryEW 
(A) 


Bry/Nal 


Bry/Hz 


Bry/Fe 11 


CO(2-0)/Nai 


HH26 IRS 


> 2 


-5.0 


1.8 


0.3 


>25 


15 


HH34 IRS 


> 5 


-8.9 


2.4 


1.7 


0.22 


11.6 


HH46 IRS 


> 1 


-4.8 


>5.1 


5.0 


0.05 


> 30 


HHIOO-IR 


6.0 ±0.5 


-7.9 


6.9 


93.0 




4.6 


IRS2 


2.9 ± 0.2 


-7.0 


19.2 


140.0 




3.5 


IRS5a 


1.0 ±0.1 


-0.6 











Notes. Columns show in order: K band veilin g, Bry equivalent w idth, Bry/Nai, Bry/H2, Bry/Fe 11 and CO(2-0)/Nai flux ratios. The same results 
for three Class I sources of the R CrA region dNisini et alj2005ah are reported for comparison (see text for details). 



hotter and internal disc region at low ionisation. We note that 
in both HH26 IRS and HH34 IRS, the FWHM velocities of the 
Na I lines are smaller than those of Bry of a factor of about two. 
This is an evidence that Bry is associated to a gas component at 
higher keplerian velocity, e.g. the accretion flow and/or ionised 
stellar/disc wind. All such findings may suggest that sources 
with energetic jets possess massive discs characterised by large 
columns of warm non-ionised gas. 

Finally, the absence of Fe 11 emission in HH26 IR points to 
the presence of a colder molecular jet in this source, as also con- 
firmed by the lower Br7/H2 ratio observed, which indicates a 
higher H2 content. 



coiTelation between amount of veiling and accretion luminosity 
has been found by N05 in the small sample of Class Is in R Cra. 

If the K band veiling is indeed due to the accretion excess 
above the stellar photosphere, then one should expect a coiTela- 
tion between the equivalent widths of accretion tracer lines and 
the veiling. This kind of correlations ha ve been indeed found fo r 
Ha and Can lines in T Tauri stars (e.g lMuzerolle et aLlll998bl) . 
It is interesting to see if such coiTelation exists also for Bry in 
the limited sample of our jet sources and CrA sources from N05. 

On the assumption that the ratio between the Bry flux and the 
K band accretion flux is roughly constant, the equivalent width 
should indeed only depend on the veiling: 



5.2. Veiling 

The lower limits we set on the K band veiling (see Tab. [8]l show 
that at least two sources have veilings greater than the ones typi- 
cally obse rved in the m ore evolved T Tauri stars, for which ~ 
1 -2 (e.g. lFolha& Emer son 1999; Johns- KruU & Val enti 2001), 
in a greement with the values measured on other Class I objects 
(e.g Doppmann et aLll2005l : iNisini et aL l l2005aLlGreene&' Lade I 
fl996r Large (rx > 2) values of the A'-band veiling are in- 
dicative of the presence of inner circumstellar dusty envelopes 
jGreene & Lada 1996) and have been suggested to be a good 
indicator of high accretion activity in the source, since ex- 
cess emission shou ld be generated by the infalling matter (e.g. 
ICalvet et al.lll997l) . DOS have measured /T-band veilings for 
Class I and flat-spectrum sources in Ophiuchus that are signifi- 
cantly higher than those of Class II in the same region and as- 
sume this as a sign of Class I stars being actually embedded ob- 
jects undergoing significant mass accretion. Finally, a positive 



EW 



"Bry rx 
K,acc ^ + fK 



(10) 



where it has been supposed that all the K band excess emission 
is due to accretion {Fi^^^c ~ Fx.acc)- From this relationship one 
should expect the line EW to increase roughly linearly for small 
rK values, reaching an asymptotic value equal to Fsry/Ffcacc for 
large values of r^. 

Despite the limited number of determinations, and the lower 
limits on rK for the jet sources, such a behaviour is approxi- 
mately recognisable in Fig.|7] where the results for the six con- 
sidered sources are plotted. In particular, the large Bry EW mea- 
sured on HH34 IRS should be noticed, which is similar to the 
value measured on HHIOO IR and in agreement with what is ex- 
pected given the large lower limit on rK- This indicates that the 
accretion luminosity in this object should indeed be high, pro- 
viding support to our estimate of an Ak value of the order of 5 
mag. 
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Fig. 7. Bry equivalent widths is plotted against the K band veil- 
ing for the three sources and the CrA sources of iNisini et alj 
(l2O05ah . 



2 4 6 8 



5.3. Accretion luminosity 

The sources here analysed have accretion luminosities that are a 
fraction between ~0.2 and 0.8 of the total luminosity. In addition 
to HH34 IRS and HH26 IRS, only three Class I have been found 
to have a Lacc/^^boi ratio > 0. 5 so far: HHlOO-IRan d IRS2 in CrA 
(N05) and YLW15 in Oph i Greene & Ladall2002h . We note that 
the fraction of Lace over Lboi in HH46 IRS is somehow uncertain, 
due to the binary nature of the object. 

Recent findings have actually shown that there seem to exist 
a large spread in the measured accretion luminosities of embed- 
ded sources. Accretion luminosities spanning from a few per- 
cent to about 50% of the total have been measured for instance 
by W04 in the Taurus-Auriga sources of their sample, although 
these results (being based on optically visible sources) would 
better apply to objects in general more evolved than our embed- 
ded tar gets. A simi lar result has been however recently found 
also bv lBeckl ( l2007h on a sample of Class I in Taurus. Finally, in 
N05 only two out of 5 investigated Class Is have bolometric lu- 
minosities dominated by accretion. D05 noted that the accretion 
luminosity derived from Bry fluxes using Muzerolle's relation 
for the sources of their sample is low (10"' < L-^cc/Lq < 10") 
even in objects with significantly high bolometric luminosities, 
which could therefore be explained only assuming a very high 
photospheric contribution. This last result could be partially 
due to an underestimate of the extinction correction to be ap- 
plied to the Bry, being the Ay value estimated from near-IR 
colours. There is nevertheless growing evidence that only a lim- 
ited number of sources normally classified as Class I actually 
have accretion-dominated luminosity. 

On the basis of the small sample we have gathered, we 
can try to tentatively define some properties characterising 
these Accretion-Dominated Young Objects (ADYOs), i.e. those 
sources with Lacc/^^boi Z 0.5. First of all, we note that all the 
identified accretion-dominated objects have rfc strictly larger 
than 2. In the D05 sample, only 6 out of 40 sources have such 
a large veiling, testifying that ADYOs might be just a fraction 
of all the Class Is. Sources with massive and dense jets detected 
very close to the exciting object have high Lacc/^^boi ratios, as we 
infer from the sample investigated in this paper, but the presence 
of a jet is not a necessary characteristic of an ADYO (e.g. CrA 
IRS 2). On the other hand, the fact that the jet sources here inves- 
tigated have SEDs typical of objects in transitions from Class 
to Class I, indicates that large accretion luminosities are associ- 
ated to less evolved sources, while a consistent number of "clas- 



sical" Class I objects may be in a phase of less intense accretion 
and have already acquired most of their final mass. 

We also note that all the ADYOs so far discovered, except 
for YLW15, have CO and Nai lines in emission, although there 
is not a direct correlation between the strength of these features 
and the accretion luminosity . The r najo rity of the Class Is inves- 
tigated in D05 as well as in iBeckl (120071) . show on the contrary 
CO and Na i features in absorption, which is a further indication 
that their samples may be biased towards sources lacking mas- 
sive accretion-heated inner discs. 

5.4. Mass accretion and ioss rates 

The actual mechanism through which accretion and ejection of 
matter take place in the protostar still remains substantially un- 
known from an observational point of view. Measurements of the 
accretion and ejection parameters are therefore of great impor- 
tance in order to provide the observational constraints needed 
to disentangle among the proposed models. A quantity of in- 
terest is the ratio Mioss/^acc, whose expected value is in gen- 
eral of the order of 0.01-0.1 in disc wind models (e.g. ,Ferreiral 
119971: ICasse & Ferr"eirall2000l) . whil e a (two-s i ded) v alue of 0.3 
is needed in the X-wind scenario of[Shuetal] ( |1994|) . 

In our case, the measured loss rates of order 10"^ - 10"^ 
Mq ' are comprised in the range of values that have b een ob- 
served for other Class I and Class II sources ( e.g. Pod io et al.l 
120061; lGullbringetaLlll998l: IWhite & Hillenbran d 2004), while 
the observed accretion rate range of 10"^ - 10"^ M0yr"' (de- 
rived from the Hi line luminosity) extends above the val- 
ues found for t he most actively accreting Class 11 jet sources 
(iMohantv et al.ll2005) . 

From these measurements we derive a mass flux that is ~0.01 
of the accretion rate in HH26 IRS and in HH34 IRS, which is in 
good agreement with values derived in T Tauri stars and with 
models, while for HH46 IRS a larger value >0. 1 is inferred (see 
TablellJ. 

Of course, the reliability of these values depend on the preci- 
sion with which Mioss and Mace are determined. These quantities 
might be affected for example by systematic errors due to the 
methods used to derive them. For instance, mass flux rates mea- 
sured from optical or near-IR line luminosities may be underes- 
timated in the presence of a dense and cold jet neutral compo- 
nent, that is not traced by the considered lines. Similarly, mass 
accretion rates estimates are affected by the intrinsic scatter of 
the Muzerolle's relation and by the stellar parameters inferred 
assuming the source on the birthline. 

Nevertheless, the observed spread could also be explained, 
at least in part, in terms of variable accretion and ejection rates 
of the sources. The accretion process (and consequently the re- 
lated ejection of matter) is in fact likely to undergo phases of 
burst and quiescence (see e.g. lLarson|[l98 4'). similarly to what 
has been suggested for other kind of accreting astronomical ob- 
jects. In this scenario Class I sources would display phases of 
enhanced accretion with a corresponding increase of the total 
luminosity. Evidences in this sense have been pr ovided for few 
Class I sources (e.g. SSV13. lLiseauetal.lll992 ') and the exis- 
tence of accretion/ejection phases is clearly observed in the jets 
themselves, where various distinct knots of emission are visible, 
often at the same distance from the source, in the red and blue 
lobes. 

We remark that variability of the accretion can produce sys- 
tematic biases when comparing the rates. This is due to the fact 
that the estimates of the mass flow can be done only on the jet 
knot closer to the protostar and therefore are derived observing 
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matter that was ejected before the accretion process ongoing at 
the moment of observations, whose signatures we detect on the 
spectrum. In this case, a better estimate would take into account, 
rather than the instantaneous accretion rate measured from the 
spectra, an accretion rate averaged over the time during which 
the matter responsible for the jet emission has been ejected. 

Unfortunately, the values of the Mioss/M-^^c ratio are proba- 
bly affected by this kind of problems in a significant number of 
cases. 

6. Conclusions 

We have investigated the accretion and ejection properties of 
three embedded sources (HH26 IRS, HH34 IRS, HH46 IRS) 
showing prominent jet-like structures. To this aim we have anal- 
ysed their medium resolution (R ~ 9000) near IR spectra ac- 
quired with VLT-ISAAC. The main results we obtained can be 
summarised as follows: 

- The bolometric luminosity and SEDs of the three sources 
have been revised on the basis of Spitzer and recent sub- 
mm observations, and theoretical radiative transfer models 
available in the literature. It turns out that the sources are 
probably very young, in a transition phase between Class 
and I. 

- The spectra of the three sources show important differences 
in their characteristics: in fact, the number and the absolute 
and relative intensity of the observed emission features (as- 
sociated both with the accretion region and the jet environ- 
ment) vary among the sources. In particular, we point out 
that there is no clear relationship between the presence of the 
jet and the spectral accretion signatures detected. Moreover, 
the spectra show no sign of absorption features, indicating 
large amount of veiling which in turn suggests the presence 
of warm dusty envelopes around the sources, heated by the 
active ongoing accretion. 

- In two of our sources (HH34 IRS and HH26 IRS) we find a 
Bry/Na i ratio much larger (by a factor of three or more) than 
that observed in other Class Is or T Tauri stars. Conversely, 
the ratio between the CO(2-0) 2.3jum overtone emission and 
the Nai 2.2Qfj.m doublet, shows the opposite trend, i.e. is 
smaller in the objects with jets. This may indicate the pres- 
ence of massive discs around the jet sources, characterised 
by large amounts of warm gas in neutral state. 

- We consistently derive Ak, L, and Lace assuming the three 
sources on the birthline and adopting the relationship be- 
tween Bry luminosity a nd accretion luminosity derived by 
iMuzerolleetaT dggSa). In the case of HH34 IRS and HH26 
IRS, accretion largely contributes to the total energy, as ex- 
pected for young sources in the main accretion phase. In 
HH46 IRS, we estimate an Lacc/Lboi ratio of the order of 
0.2, but the real value largely depends on how the bolomet- 
ric luminosity is shared among the two binary components. 
Taking into account /T-band veiling values rfc greater than 
the lower limits inferred from the spectra would lead in gen- 
eral to higher A /f. Lace and Mace- 

- Mass accretion and loss rates span (including errors) in the 
range lO^'^-lO"'' M0yr"'. The derived MhsJMacc ratio is 
~0.01-0.03 for HH26 IRS and HH34 IRS, and >0.1 for 
HH46 IRS. These numbers are in the range of values usually 
predicted by models and found in the more active classical T 
Tauri stars. 

- Comparing the results found in this work with other spec- 
troscopic studies recently performed on Class I sources. 



we conclude that the number of Class Is actually hav- 
ing accretion-dominated luminosities (Accretion-Dominated 
Young Objects, ADYOs) could be limited. From the proper- 
ties inferred in a small sample of objects we can tentatively 
define some criteria to characterise such sources: ADYOs 
have all /T-band veiling larger than 2 and in the majority of 
the cases present (in addition to H i) IR features of CO and 
Na I in emission, although these latter do not directly corre- 
late with the accretion luminosity. Class Is with massive jets 
have high Laee/Lboi ratios but not all the identified ADYOs 
present a jet. The SEDs of our small sample of three objects, 
suggest that accretion-dominated sources could be in an evo- 
lutionary phase in transition between Class and I. 
Of course, studies of the kind presented here but carried 
out on larger samples of possible candidates should be per- 
formed in order to test and refine the tentative criteria that 
we have just mentioned. 
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